function [] = fig_sec10_5_csmaBackoffOptimizationDual()

close all;

ue = load ('data/csmaBackoffOptimizationDual_ue.txt');
re = load ('data/csmaBackoffOptimizationDual_re.txt');
netUtilityWithoutBeta = load ('data/csmaBackoffOptimizationDual_netUtilityWithoutBeta.txt');
netUtilityWithBeta = load ('data/csmaBackoffOptimizationDual_netUtilityWithBeta.txt');
objFunc = load ('data/csmaBackoffOptimizationDual_objFunc.txt');

jom_ue = load ('data/csmaBackoffOptimizationDual_jom_ue.txt');
jom_re = load ('data/csmaBackoffOptimizationDual_jom_re.txt');
jom_netUtilityWithoutBeta = load ('data/csmaBackoffOptimizationDual_jom_netUtilityWithoutBeta.txt');
jom_netUtilityWithBeta = load ('data/csmaBackoffOptimizationDual_jom_netUtilityWithBeta.txt');
jom_objFunc = load ('data/csmaBackoffOptimizationDual_jom_objFunc.txt');

% multiplot ('Time', 'Objective function' , objFunc-netUtilityWithBeta+netUtilityWithoutBeta , jom_objFunc-jom_netUtilityWithBeta+jom_netUtilityWithoutBeta , [-100 -30] , '-' , 'r_e' , re , jom_re , [0 150] , '-' , 'u_e' , ue , jom_ue , [0 0.5] , '-');
% multiplot ('Time', 'Objective function' , objFunc , jom_objFunc , [-100 -30] , '-' , 'r_e' , re , jom_re , [0 150] , '-' , 'u_e' , ue , jom_ue , [0 0.5] , '-');
% multiplot ('Time', 'Objective function' , netUtilityWithoutBeta , jom_netUtilityWithoutBeta , [-100 -30] , '-' , 'r_e' , re , jom_re , [0 150] , '-' , 'u_e' , ue , jom_ue , [0 0.5] , '-');
multiplot ('Time', 'Objective function' , netUtilityWithoutBeta , jom_netUtilityWithoutBeta , [-100 -30] , '-' , 'r_e' , re , jom_re , [0 1500] , '-' , 'u_e' , ue , jom_ue , [0 0.5] , '-');
% multiplot ('Time', 'Objective function' , netUtilityWithBeta , jom_netUtilityWithBeta , [-100 -50] , '-' , 'r_e' , re , jom_re , [0 150] , '-' , 'u_e' , ue , jom_ue , [0 0.5] , '-');


jom_objFunc
grad_objFunc = objFunc (end)

jom_netUtilityWithoutBeta
grad_netUtilityWithoutBeta = netUtilityWithoutBeta (end)

jom_netUtilityWithBeta
grad_netUtilityWithBeta = netUtilityWithBeta (end)

max(ue(end,2:end))
max(jom_ue)

